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Q\ ■ We introduce a variable useful for measuring masses of particles pair produced at 

(— I . hadron colliders, where each particle decays to one particle that is directly observable 

Qh! and another particle whose existence can only be inferred from missing transverse 

Q^' momenta. This variable is closely related to the transverse mass variable commonly 

^ • used for measuring the W mass at hadron colliders, and like the transverse mass our 

variable extracts masses in a reasonably model independent way. Without considering 

either backgrounds or measurement errors we consider how our variable would perform 

^ ■ measuring the mass of selectrons in a mSUGRA SUSY model at the LHC. 
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Before 1983, e+e" colliders dominated the search for new particles, their clean final states 
making identification of new particles far easier than at hadron machines. However in 1983 
this changed with the observation of the W and Z bosons by the UAl||l], |^ and UA2[0, ^ 
collaborations at the CERN pp collider. Despite the more difficult conditions where the 
initial state protons and anti-protons gave rise to an underlying event which complicates 
measurement of the final state, the higher energies available to a pp machine meant they 
discovered the W and Z bosons before other machines had enough energy to produce them. 
This trend continued in 1994/5 when CDF and D0 discovered the t quark at the Tevatron 
pp collider 0. To date the only direct experimental evidence that we have for the t quark 
is from the Tevatron. 

Clearly, the greater energies available to pp machines gives them an advantage in discov- 
ering fundamental particles, and we can hope that this trend will be continued at the LHC. 
Upon discovery of a new particle the first question that is asked is what is the mass of this 
particle. There are two fundamentally different ways of measuring the mass of a particle. 
Usually one tries to extract the mass of the particle directly from the observed events. For 
example UAl and UA2 detected the Z boson in its decay to lepton pairs e~^e~ or jX^n' . For 
each event that contains a pair of leptons one can construct the mass of the particle that 
produced the leptons as, 

w? = {pi++pi-f (1) 

and on an event by event basis obtain a direct estimate the mass of the particle. Alternatively 
one could make an indirect model dependent measure of the Z mass by measuring the 
"Z event" cross-section, and using a model that predicts the Z cross-section in terms of the 
Z mass, perform a fit to the measured cross-section. For measuring the Z mass the direct 
measurement is superior, if only for its lack of model dependence, however when measuring 
the mass of the W boson the case is less clear. At the CERN pp collider the W boson 
was detected in its decay to a charged lepton and a neutrino, however the neutrino goes 
unobserved and so only gives rise to missing momentum. This means that one can't directly 
observe the W mass from the lepton and neutrino momenta. UAl [Q formed the transverse 
mass variable, 

ml = 2{E'j.$T-VT-i>T), (2) 

where ^1^ = fi^. This variable has the property that, 

m]. < m^ (3) 

with equality possible for events where the lepton and neutrino are produced with the same 
rapidity. Thus although one can't obtain the W mass from a single event, one can obtain 
a lower limit on the W mass. In addition if one obtains the lower mass bound from many 
events, this can approach the W mass, and so one can extract the W mass in a model 
independent way. UAl |jl| along with UA2 f^ also performed a model depenendent fit to 
the pt spectrum for the lepton, extracting the W mass from that model. 

In practice, though, measuring the W mass with ttit does have some small model depen- 
dence, as the precise fraction of events which occur with rriT close to rriw is still dependent 




Figure 1: Diagram of the generic process that we consider. A hadronic colhsion that leads to 
a pair of particles being produced, which each decay into one particle that is observed with 
momenta pi and p2 respectively; and one particle (shown as a wavy lines) that is not directly 
detected, and whose presence can only be inferred from the missing transverse momentum, 

on the physics processes which produce the W boson, and how the W boson decays. In 
addition, the missing transverse momentum is poorly measured experimentally compared 
with prii), so the theoretical model dependence of the measurement of niw from the pr(0 
spectrum is balanced by the experimental error on extracting mw from the edge of the rriT 
spectrum. 

In this paper we wish to introduce a variable which measures particle masses, which like 
transverse mass has little dependence on exactly how such massive particles are produced. 
The variable is used for the generic process shown in figure |I|, where a hadronic collision pair 
produces a massive particle whose dominant decay is into one observed and one unobserved 
particle. This unobserved particle can only be detected from the missing momentum that it 
carries away, and that the massive particle is pair produced means that we can only measure 
the missing momentum of the pair of invisible particles. Although this may sound like an 
unusual process to look for new particles, it naturally occurs in any theory where there is an 
(approximately) conserved charge, and the lightest particle with that charge is only weakly 
interacting. Two examples of where such a situation can occur are SUSY models and a 
4th lepton generation. In R-parity conserving SUSY models, sparticles are pair produced, 
and cascade decay to the lightest sparticle, which must be stable and is expected to not be 
directly detectable. Slepton production and decay can often follow this route: 

pp^x + ip], ^x + i^rxiXv (4) 

In such an event the pair of lightest SUSY particles, x? go unobserved and only leave their 
signature as missing transverse momentum. 

For a 4th generation lepton the charged lepton would be pair produced in a Drell-Yan 
type process, decaying to a neutrino and a W boson, 

pp^X + Itll ^ X + uuW+ui,W- (5) 



the (probably massive) ui^ going unobserved, while the W bosons could be detected in their 
decays to either lu or to jets. 

We now look specifically at the process given in equation ^, although the variable which 
we define would work identically in any process where a particle is pair produced and decays 
to one visible and one invisible particle. 

The variable that we wish to introduce is closely related to my, however the standard 
definition of my, given in equation ||, assumes that the unobserved particle is massless, so 
we return to the derivation of this variable. For the decay, 

l^lx (6) 

for arbitrary momenta we can write, 

m? = mf + m| + 2{EtiEtx cosh(Ar7) - prz ■ Ptx) (7) 

where Et = \/'Pt + ^ ^"^^ '^V is the difference in rapidity, '>] = \ ln[(_E + pz)/{E — pz)], 
between between the / and x. 

Now as coshr] > 1 we have, 

m? > m^(pT;, Ptx) = rnf + ^l + '^{EtiEtx - Vti ■ Ptx)- (8) 

This gives a version of transverse mass valid for arbitrary masses, with equality when the / 
and X ^'^^ produced with the same rapidity. Notice that Eti and Etx depend on m^ and m| 
respectively. 

The transverse mass can't be formed directly from the process in equation (H), as both 
the neutralinos give rise to missing momentum, however we can experimentally measure the 
sum of their transverse momenta as the missing transverse momenta in the event, 

f^T = PTxa+PTx6- (9) 

If PTxa ^"^^ Ptx6 were obtainable, then one could form two transverse masses, and using the 
relationship (|D obtain, 

m? > max{m|(pT;-,PTxa),"^T(Pr;+,PTxJ} (10) 

However, not knowing the form of the splitting (^, the best we can say is that: 

max {m|(pr«- , jii), m|(pTi+, 5^2)} (H) 






With the minimization over all possible 2-momenta, 1^1,2, such that their sum gives the 
observed missing transverse momentum, ]^t- This is the variable, called Mt2, that we wish 
to introduce. This bound we can obtain directly from experimentally measured parameters. 
Although not totally transparent, for particular momenta, Mt2 can be equal to mf the 
requirement being that for both slepton decays the lepton and neutralino are produced at 



the same rapidity (although the sleptons themselves can be at differing rapidities), and in 
addition, 

J^Tl- ^Txa^ ^^Tl+ ^Txb^ 

We have not managed to derive a general analytic expression for the minimization over 
splittings of T^T\ largely because an experimental measurement of piy only measures the 
missing transverse momentum, and neither the missing energy nor the missing longitudinal 
momenta. This means that ^t is not a 4 vector, which means that Mj'2 can not be calculated 
in a manifestly Lorentz invariant manner. The complication that this introduces to the 
minimization is enough to making an analytic solution non trivial. However if we take one 
of the parameters that we minimize over to be rn^ = j^^ {i.e. that ^t = \/fTi'^ + i>x) 5 ^^nd 
in addition note that longitudinal momenta play no part in the definition of Mx2, then the 
remaining minimization becomes Lorentz invariant, and this means that it can be solved 
analytically. This leaves one minimization over m^ that we performed numerically. The 
form of this minimization is not particularly illuminating, and so we do not give it here; 
however a computer code for evaluating M^2 is available from the authors. 

Of course, the variable Mt2 is only good at extracting particle masses from processes 
having many events close to the maximum allowed value, and this depends on the physics 
of the process being measured. Hence, a priori, we can not say that Mt2 is useful for 
measuring particle masses in all processes. However we expect typical physics processes to 
have reasonable numbers of events close to the maximum value allowed for Mjn2 and hence 
expect it to be a useful variable. To illustrate the variable in use, we consider a model 
which allows the process shown in equation (^, at the expected LHC centre of mass energy 
^/s = 14TeV. Our SUSY model is the fifth minimal supergravity model (mSUGRA) ^ 
point selected by the LHC Committee in 1996 for detailed study by the ATLAS and CMS 
collaborations 0. This model is characterized by, 

{tan/5 = 2.1, mi/2 = 300 GeV, mo = 100 GeV, Aq = 300 GeV, fi > 0}. (13) 

We generate events using the SUSY version of the Herwig generator [§ . For our purposes 
the only important features of the model that we use to produce our events are, 

rrii^ = 157.1 GeV m^o = 121.5 GeV. (14) 

The angular distributions for the sleptons are given by the Drell-Yan process that produces 
spin-0 particles, while the decay of the sleptons is isotropic due again to their being spinless. 
Herwig serves to build up a reasonably realistic underlying event, i.e. that pi^ 7^ Pti+ +Pti-- 
Defining the missing momenta as, 

i>T = PTxa + PTx„ (15) 

we generate 1105 events, which corresponds to an integrated luminosity of about 30 fb^^, 
which the LHC should collect in approximately one year of low luminosity running. 
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Figure 2: Mt2 distribution for the process pp 
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X + 1+1- XiXi at the LHC. With 



mj = 157.1 GeV and m^ = 121.5 GeV, assuming the actual value for m^ when calculating 
Mj'2- The data with error bars are 1105 events, that the LHC would collect in approximately 
1 year of running at low luminosity, i.e. £ ~ 30 fb~^. The histogram represents C ^ 500 fb~^ 
to show the shape of the distribution that would be obtained with huge statistics, with the 
normalization modified to be the same as £ ~ 30 fb~^. 
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Figure 3: Values of mj that would be obtained from that largest Mt2 value observed, where 
differing values of m^ are used in the calculation of Mt2, for the 1105 events shown in figure 
1^. All events have been generated with mi = 157.1 GeV and m^ = 121.5 GeV. 



In figure ^ we show the Mj'2 distribution for this modeL Quite clearly the maximum 
Mt2 value allowed corresponds well to the mass of the selectron, with the Mt2 distribution 
tending smoothly to zero at that point. While such an edge is not as easy to measure as a 
vertical drop, we still expect that such an edge would be detectable experimentally, and hence 
a viable means to measure the selectron mass. The lower bound on Mt2 is rn^^ however this 
does not give a means to measure the LSP mass, the lower bound on Mt2 follows directly 
from equation (^, and this is a parameter that goes into the calculation of Mt2- Indeed if 
the value of m-^ is unknown in addition to mi one can only obtain a relationship between 
the two. In figure ^ we show the value of mi that would have been obtained from the same 
events shown in figure |^, using differing input values for m^ in the calculation of Mt2- One 
can see that slepton mass that one would extract from the events behaves approximatly as, 

m^i ~ m,^ + constant, (16) 

and so any uncertainty in the LSP mass will be directly reflected as an error in the extracted 
selectron mass. It should also be noted that the edge of the Mt2 distribution becomes hard 
to fit as the input x mass deviates from the physical x mass as the Mt2 distributions develops 
a tail at the largest Mt2 values. 

To conclude, in this paper we have introduced a new variable for measuring masses 
of particles produced at hadron colliders, where the longitudinal momentum of the hard 
scattering is typically unmeasured. It may be used when particles are pair produced, with 
each decaying to one particle that is directly observed and one particle that is not directly 
observed. This variable is analogous to the transverse mass variable, m^, commonly used 
for measuring the W mass in its decay to Iv at pp colliders, except that it works where 
the particle being measured is pair produced and where the unseen particle is massive. We 
expect that the masses extracted using our variable, like those from the transverse mass, 
will be largely independent of the physics processes which produce the particles, and hence 
give a viable means of extracting masses in a model independent way. As an illustration of 
this variable in action we consider measuring the mass of selectrons in SUSY at the LHC 
The results look promising however as we have not considered the effects of any background 
processes or experimental mis-measurement errors further study is required[0. 
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